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Abstract 

The gamma-ray burst (GRB) 050904 at z = 6.3 provides the first opportunity of probing the intergalactic 
medium (IGM) by GRBs at the epoch of the reionization. Here we present a spectral modeling analysis 
of the optical afterglow spectrum taken by the Subaru Telescope, aiming to constrain the reionization 
history. The spectrum shows a clear damping wing at wavelengths redward of the Lyman break, and the 
wing shape can be fit either by a damped Lya system with a column density of log(iVHi/cm“^) ~ 21.6 at a 
redshift close to the detected metal absorption lines (Zmetai = 6.295), or by almost neutral IGM extending 
to a slightly higher redshift of Zigm,u ~ 6.36. In the latter case, the difference from Zmetai may be explained 
by acceleration of metal absorbing shells by the activities of the GRB or its progenitor. However, we 
exclude this possibility by using the light transmission feature around the Ly/3 resonance, leading to a firm 
upper limit of Zigm,u < 6.314. We then show an evidence that the IGM was largely ionized already at 
z = 6.3, with the best-fit neutral fraction of IGM, a;Hi = 0.00, and upper limits of xhi < 0.17 and 0.60 at 68 
and 95% C.L., respectively. This is the hrst direct and quantitative upper limit on xhi at z > 6. Various 
systematic uncertainties are examined, but none of them appears large enough to change this conclusion. 

To get further information on the reionization, it is important to increase the sample size of z ^ 6 GRBs, in 
order to find GRBs with low column densities (logiVni ^ 20) within their host galaxies, and for statistical 
studies of Lya line emission from host galaxies. 

Key words: cosmology: early universe — gamma rays: observations — gamma rays: theory — 
galaxies: intergalactic medium 


1. Introduction 

Although more than 30 years have passed since the dis¬ 
covery of gamma-ray bursts (GRBs) (Klebesadel, Strong, 
& Olson 1973, see, e.g., Meszaros 2002; Piran 2004 for 
recent reviews), it was rather recent that GRBs became 
widely recognized as a unique tool of cosmological studies 
and exploring the early universe. It was just before GRBs 
were proven to have the cosmological origin (Metzger et 
al. 1997) that the first attempt to study the cosmic star 
formation history by using GRBs was made (Totani 1997; 
Wijers et al. 1998). The potential use of GRBs as a probe 
of the reionization history of the intergalactic medium 
(IGM) was also pointed out (Miralda-Escude 1998), but 
at that time, it was thought that GRBs can probe the 
universe at most up to modest redshifts of z ~ 3 by instru¬ 
ments available at that time or in the near future. Such 
notion was, however, soon discarded by the discoveries of 
extremely luminous GRBs like GRB 971214 (Kulkarni et 
al. 1998) and 990123 (Kulkarni et al. 1999), which could 

* Based on data collected at the Subaru Telescope, which is oper¬ 
ated by the National Astronomical Observatory of Japan. 


be detected even at redshifts beyond 10. Then it did not 
take long time before astrophysicists started to discuss 
GRBs as a promising lighthouse to study the extremely 
high redshift universe, potentially giving cosmologically 
important information including the population HI star 
formation and the reionization history (Lamb & Reichart 
2000; Ciardi & Loeb 2000). 

During 2000-2004, satellites such as the BeppoSAX 
and the HETE-2 continued to discover more and more 
GRBs, and there was important progress including the 
establishment of the firm connection between long dura¬ 
tion GRBs and energetic supernovae (Stanek et al. 2003; 
Hjorth et al. 2003a). However, the distance of GRBs 
did not extend to very high redshift, with the highest 
record of z = 4.5 (Andersen et al. 2000). The launch 
of the Swift satellite in 2004 then allowed the GRB com¬ 
munity to search for fainter and more distant GRBs with 
improved detection rate. The GRB 050904 was discov¬ 
ered by the Swift on 2005 September 4 at 01:51:44 UT 
(Gusumano et al. 2005), and follow-up photometric ob¬ 
servations of the afterglow found a strong spectral break 
between optical and near-infrared (NIR) bands, indicat¬ 
ing a very high redshift of z ~ 6 (Haislip et al. 2005; Price 
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et al. 2005; Tagliaferri et al. 2005). This suggestion was 
confirmed by the subsequent spectroscopic observation by 
the Subaru Telescope, which found metal absorption lines 
at z = 6.295 and the corresponding Lyman break and red 
damping wing (Kawai et al. 2005, hereafter Paper I). This 
opens a new era of GRB observations at redshifts that are 
close to the cosmic reionization and comparable to those 
of the most distant galaxies (Taniguchi et al. 2005) and 
quasars (White et al. 2003; Fan et al. 2006). Here we re¬ 
port a detailed analysis and interpretation of the optical 
afterglow spectrum of GRB 050904 presented in Paper I, 
to derive implications for the reionization. 

The famous Gunn-Peterson (GP) test tells us that the 
IGM is highly ionized at z ^ 5 (Gunn & Peterson 1965), 
while the observations of the cosmic microwave back¬ 
ground radiation (GMB) indicates that the universe be¬ 
came neutral at the recombination epoch of z ~ 1100 
(Spergel et al. 2003). The reionization of the IGM is 
believed to have occurred during z ~ 6-20 by the first 
stars and/or quasars, and the precise epoch and nature of 
the reionization is one of the central topics in the modern 
cosmology (see Loeb & Barkana 2001; Barkana & Loeb 
2001; Miralda-Escude 2003; Haiman 2004 for recent re¬ 
views) . The dramatic increase of the optical depth of the 
Lya forest with increasing redshift at z ^ 5.2 (Djorgovski 
et al. 2001) and the subsequent discovery of broad and 
black troughs of Lya absorption (the GP troughs) in the 
spectra of z > 6 quasars (Becker et al. 2001; White et al. 
2003; Fan et al. 2003, 2006) indicate that we are begin¬ 
ning to probe the epoch of reionization (Fan et al. 2002; 
Gen & McDonald 2002; but see also Songaila & Gowie 
2002). On the other hand, the polarization observation of 
the GMB by the Wilkinson Microwave Anisotropy Probe 
(WMAP) indicates a much higher redshift of reionization, 
z = 17±5 (Kogut et al. 2003). Some theorists have argued 
that the hydrogen in the IGM could have been reionized 
twice (Wyithe & Loeb 2003; Cen 2003). 

Because the cross section of the Lya resonance absorp¬ 
tion is so large, the light blueward of the Lya wavelength 
at the source is completely attenuated if the IGM neu¬ 
tral fraction ihi = nm/nn is larger than ~ 10“^, and 
hence the Lya trough of z ~ 6 quasars gives a constraint of 
only a;Hi ^ 10“^. The cross section becomes much smaller 
for longer wavelength photons than the Lya resonance, 
and the spectral shape of the red damping wing of the 
GP trough can potentially be used to measure xhi more 
precisely (Miralda-Escude 1998). However, applying this 
method to quasars is problematic because of the uncer¬ 
tainties in the original unabsorbed quasar spectra and the 
proximity effect (Bajtlik, Duncan, & Ostriker 1988), i.e., 
the ionization of surrounding IGM by strong ionizing flux 
from quasars (Gen & Haiman 2000; Madau & Rees 2000). 
Therefore, though some authors suggested xui ^0.1 using 
z ~ 6 quasar spectra (Mesinger & Haiman 2004; Wyithe & 
Loeb 2004; Wyithe, Loeb, & Garilli 2005), these estimates 
are generally model dependent. 

The Lya line emission is seriously attenuated if it is em¬ 
bedded in the neutral IGM, and hence the Lya line emis- 
sivity of galaxies at z ^ 6 is another probe of the reioniza¬ 
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tion. Therefore the detection of Lya emission from many 
galaxies at z ^ 6 (Hu et al. 2002; Kodaira et al. 2003; 
Taniguchi et al 2005) may indicate that the universe had 
already been largely ionized at that time (Malhotra & 
Rhoads 2004; Stern et al. 2005; Haiman & Gen 2005). 
However it should not be naively interpreted as imply¬ 
ing small xhi, since these Lya emitters (LAEs) are se¬ 
lected by strong Lya emission and hence they may be 
biased to those in ionized bubbles created by themselves 
or clusters of undetected sources (Haiman 2002; Wyithe & 
Loeb 2005). On the other hand, the Lyman-break galax¬ 
ies (LBGs) selected by broad-band colors are free from 
the selection bias about Lya emission, but the statistical 
nature of Lya emission from LBGs at z ^ 5 is not yet 
well understood (Ando et al. 2004; Bouwens et al. 2004; 
Dickinson et al. 2004; Giavalisco et al. 2004; Stanway 
et al. 2004), compared with those at z ~ 3 (Shapley 
et al. 2003). Furthermore, the Lya line emission from 
lower-redshift starburst galaxies is often redshifted with 
respect to the systemic velocity of galaxies (e.g., Pettini 
et al. 2000), and such a relative redshift will lead to a 
higher detectability of LAEs at z ^ 6, indicating a pos¬ 
sible systematic uncertainty in the reionization study by 
Lya emission. 

GRBs have a few advantages as a probe of the cos¬ 
mic reionization, compared with quasars or LAEs/LBGs. 
GRB afterglows are much brighter than LAEs/LBGs and 
comparable to or even brighter than quasars if they are 
observed quickly enough after the explosion. The Lya or 
ultraviolet luminosity of host galaxies is irrelevant to the 
detectability of GRBs, and hence GRBs can probe less 
biased regions in the early universe, while quasars and 
bright LAEs/LBGs are likely biased to regions of rapid 
structure formation with strong clustering. In most cases 
it is expected that the IGM ionization state around GRB 
host galaxies had not yet been altered by strong ionizing 
flux from quasars. Finally, the spectrum of GRB after¬ 
glows has a much simpler power-law shape than compli¬ 
cated lines and continuum of quasars and LAEs/LBGs, 
and hence model uncertainty can greatly be reduced. 
Especially, a detailed fitting analysis of the damping wing 
of the GP trough in an afterglow spectrum may lead to a 
precise determination of xhi. 

In this paper we take these advantages for the first time 
and derive the first implications for the reionization from 
GRBs by using the spectrum of GRB 050904. In §2, we 
briefly describe the overall features of the observed spec¬ 
trum. Formulations of the model fitting will be given in 
§3. Before the fitting, theoretically possible ranges of some 
model parameters will be defined (§4). Following the de¬ 
scription of the fitting procedure (§5.1), the fitting results 
to the red damping wing will be given in §5.2 and §5.3. 
The additional constraint from the Ly/3 feature is dis¬ 
cussed in §5.4. We then derive a quantitative constraint 
on the IGM neutral fraction, taking into account all con¬ 
straints and identified uncertainties in §6. In §7 we discuss 
the prospects of the future GRB data as the reionization 
probe, based on the lessons from this first data set; §7.1 
is for the neutral hydrogen column density within host 
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galaxies, and §7.2 is for the Lya line emission from host 
galaxies. A summary and conclusions will be presented in 
§ 8 - 

Throughout this paper, we use the WMAP values of 
the cosmological parameters in the flat universe: Hq = 
71 km s“^Mpc“^, = 0.044, and VIm = 0.27 (Spergel et 

al. 2003), and the primordial helium mass fraction in the 
total cosmic baryon, Yp = 0.25 (Kawasaki, Kohri, & Moroi 
2005). 

2. The Afterglow Spectrum of GRB 050904 

Figure 1 presents the spectrum of GRB 050904 at 
wavelengths longer than the Lyman break, taken by the 
FOGAS instrument (Kashikawa et al. 2002) mounted on 
the Subaru Telescope (lye et al. 2004). Figure 2 presents 
the close-ups of the Lya and Ly/3 resonance regions. See 
Paper I for the spectrum in the full wavelength range, 
the details of the observation, and basic results including 
the list of line detections. The pixel scale is 2.67 A/pix 
and the wavelength resolution of the spectrum is 8.5 A 
(FWHM) at Aobs 9000A. The total exposure time is 4.0 
hours and the mid time of the exposure was September 7, 
12:05 UT (3.4 days after the burst). 

A clear damping wing is seen only redward of Aobs ~ 
8900 A, corresponding to the Lya wavelength Aq, = cjva = 
1215.67A at the redshift of the identified metal absorption 
lines, Zmetai = 6.295. All the metal lines were identified 
at the same redshift within the error of ± 0.002, except 
for an intervening absorption line system at z = 4.840. 
We searched any time variability of the damping wing 
shape and the absorption line strength by separating the 
spectrum into the first 1.5 hrs, the mid 1 hr, and the 
final 1.5 hrs, but statistically significant variability was 
not found. 

At wavelengths redward of the Lya, the spectral en¬ 
ergy distribution (SED) has almost no feature except for 
the damping wing and the metal absorption lines, as ex¬ 
pected from the power-law SED of GRB afterglows. There 
are some apparent emission-like features at Aobs ~ 8820- 
8920A, which may be Lya emission. However, the data 
of this wavelength range is relatively noisy and we cannot 
confirm these signals to be real. (See the error and the 
CGD image of the spectrum shown in Fig. 2.) We thus 
set only a conservative upper limit on the possible Lya 
emission (§7.2). 

In wavelengths blueward of the Lya, the spectrum re¬ 
markably resembles those of quasars at similar redshifts. 
The deep, black Lya trough can be seen in Aobs ~ 8500- 
8900 A where the flux level is consistent with zero. On the 
other hand, transmission of the afterglow light can be seen 
in ~ 7500-8500 A (Paper I), implying that the universe 
became more transparent against the Lya absorption at 
z ^ 6. The transmission, however, again disappears at 
Aobs ^ 7500 A corresponding to the expected trough by 
the Ly/3 absorption. 


3. Spectral Modeling 

The clear red damping wing tells us the existence of a 
large amount of neutral hydrogen along the line of sight to 
the GRB. There are two possibilities for the location of the 
neutral hydrogen. One is a damped Lya system (DLA) 
associated with the host galaxy, as often seen in the opti¬ 
cal spectra of high redshift GRB afterglows (Jensen et al. 
2001; Fynbo et al. 2001; Mpller et al. 2002; Hjorth et al. 
2003b; Vreeswijk et al. 2004; Ghen et al. 2005; Starling 
et al. 2005; Vreeswijk et al. 2005; Watson et al. 2005b; 
Berger et al. 2005a). The other is the absorption by neu¬ 
tral hydrogen in the IGM before the completion of the 
reionization, as discussed in §1. Here we try to fit these 
two components to the observed shape of the damping 
wing, and constrain the parameters of the neutral hydro¬ 
gen column density of the DLA, N^i (in cm“^), and the 
neutral hydrogen fraction in the IGM, xhi- 

We assume the original afterglow spectrum before ab¬ 
sorptions as iG oc with the intrinsic power index of /3o • 
There are a few reports of jSohs in the NIR bands based on 
color observations: /3obs = —10.6 hrs (Haislip 
et al. 2005), /3obs = —0.3 ±0.6 at 0.51 days (Price et al. 
2005), and ,5obs = — 1-25±0.25 at 1.155 days (Tagliaferri et 
al. 2005). The NIR colors do not show any indication of 
extinction at the host galaxy (Tagliaferri et al. 2005). The 
X-ray data taken by the Swift shows a large soft X-ray ab¬ 
sorption corresponding to log Ahi = 22.86 in the first ~250 
s from the trigger, but it disappears from the data after 
~500 s with an upper bound of logAni < 21.57 (Boer et al. 
2005), corresponding to a visual extinction of Ay < 2.1 by 
the standard correlation for the solar abundance (Predehl 
& Schmitt 1995)^. Later in this paper we will find that 
logAni ^ 21.6 from the fitting to the observed damping 
wing shape, and metallicity of ^ ~ O. 1 (Ahi/ 1 O^^ ®)“^Z 0 
from the SlI column density. By using the relation inferred 
for the Small Magellanic Cloud (SMC) having a similar 
metallicity (Bouchet et al. 1985), the Hi column density 
translates to Ay < 0.3. In fact, the observed reddening 
of optical GRB afterglows is generally even smaller than 
those expected from the corresponding X-ray absorption 
or column densities of Hi or metals (Galama & Wijers 
2001; Savaglio & Fall 2004; Watson et al. 2005b; De 
Pasquale et al. 2005), possibly due to sublimation of dust 
grains by strong GRB radiation (Waxman & Draine 2000; 
Fruchter et al. 2001). Therefore we assume Ay = 0 at the 
host galaxy in our baseline model. 

According to the standard synchrotron shock model of 
GRB afterglows (Sari, Piran, & Narayan 1998), the ob¬ 
served index is rather soft corresponding to the electron 
spectral index oi p = 2/3o — 1 ~ —3.5, indicating that the 
frequency of the cooling break in the afterglow SED had 

^ The column densities derived by Boer et al. (2005) are also 
assuming the solar abundance (M. Boer, a private communica¬ 
tion). Since the soft X-ray absorption is dominated by metals, 
this conversion is not sensitive to a hydrogen-to-metal conver¬ 
sion (e.g., Watson et al. 2005b). On the other hand, the detec¬ 
tion of X-ray absorption for GRB 050904 may not be significant 
(see Watson et al. 2005a). 
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already passed through the NIR bands to longer wave¬ 
bands. Therefore we do not expect a spectral change be¬ 
tween the epochs of the photometric observations and our 
spectroscopic observation. We then take /3o = —1.25 in 
our baseline model. However, both the change of Pq and 
the effect of possible reddening by dust in the host galaxy 
will be discussed later. 

The optical depth of the DLA is calculated in the stan¬ 
dard manner: rDLA(Aobs) = 1 Vhi crQKbs(l + ^DLA)], where 
r'obs = c/Aobs is the observed frequency and zdla is the 
redshift of the DLA. The damping wing shape is much 
wider than the Doppler broadening by a reasonable veloc¬ 
ity dispersion (~ 3 A for 100 km/s), and hence we ignore 
it. We use the exact formula of the Lya absorption cross 
section for a restframe frequency v. 

8^ 47r2(i/-z.„)2 + A2(i./i/„)6/4 ’ 

where fa = 0.4162 and Aq = 3{gu/gi)~^ fa^ci,a are the ab¬ 
sorption oscillator strength and the damping constant of 
the Lya resonance (Peebles 1993; Madau & Rees 2000), 
respectively, and gu and gi [gu/gi = 3 for Lya) are the 
statistical weights for the upper and lower levels, respec¬ 
tively (Cox 2000). Here, the classical damping constant is 
Aci,a = ( 87 r 2 e^)/( 3 mecA^) = 1.503 x 10® s“^, where e and 
TOe are the electron charge and mass, respectively. 

For the optical depth of the red damping wing by IGM 
absorption, we use the formula given by Miralda-Escude 
(1998): 


XhI Ra TGp(^host) 

■PlGMiAobsj = - 

TT 




■ 


1 “t“ 2 :host / 

1 + ^:iGM,u \ _ j 1 + ^IGM,l' 


1 T - 2 ^obs 


1 “t” 2:obs 


where (1 + Zobs) = Aobs/Aa, i?a = AQAQ/(47rc) = 2.02x 10~® 
and Zhost is the redshift of the GRB host galaxy. In this 
formula, the IGM is assumed to be uniformly distributed 
in a redshift range from zigm,i to zigm.u- Since the damp¬ 
ing wing is by definition optically thin, the dumpiness of 
the IGM does not have significant effect. We can get a 
reasonable constraint on the mass-weighted optical detph 
by assuming a uniform distribution, which is another ad¬ 
vantage of using the damping wing. The function I{x) is 
given by: 


^9/2 Q q 

J(^) = ^ + + 3x®/® + 

l-x 7 5 

9 


( 3 ) 


which is almost always a sufficiently good approximation 
under the condition of (Zobs —-Zigm.u) ^-Ra(l+ -Zobs)- The 
Gunn-Peterson optical depth is given as: 


,, ^faAcl,aXlPc^B{l-Yp)^^ 

tgp(^) =--775-(1 ■ 


STTrupHonli^ 


.)3/2 


= 3.88 X 10® 


3/2 


( 4 ) 

( 5 ) 
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where nip is the proton mass and pc = 3Hq /(SttG) is the 
standard critical density of the universe at the present 
time. The Icr uncertainty of this numerical factor coming 
from the estimated errors of the cosmological parameters 
(Spergel et al. 2003; Kawasaki, Kohri, & Moroi 2005) is 
less than 8 %. It should be noted that, even though the 
dependence on Zhost appears explicitly in eq. ( 2 ) by the 
normalization of tgp at the host redshift, it actually can¬ 
cels out with the factor (H-Zhost)”®'^^ in eq. ( 2 ) and hence 
■PiGM depends only on zigm,i and zigm.u- Therefore Zhost is 
not included in our fitting parameters. We set zigm,i = 6 , 
and dependence on this parameter will be discussed later. 
The possible range of zigm u will be discussed in detail in 

§ 4 . 

Although the neutral fraction xm cannot exceed the 
unity, what is relevant for the IGM absorption is the ab¬ 
solute density of neutral hydrogen rather than the neu¬ 
tral fraction, and the total cosmic hydrogen density de¬ 
pends on the cosmological parameters. Therefore, we al¬ 
low xhi > 1 in this paper, as a parameter meaning the 
neutral hydrogen density normalized by the total hydro¬ 
gen density estimated from the standard values of the cos¬ 
mological parameters. 

4. Possible Ranges of zdla and 2 igm,u 

Now we have the primary fitting parameters of iVui, 
-2DLA, a:Hi, and Zigm,u- Though zdla and zigm.u can be 
treated as free parameters in the fitting, metal absorption 
lines at Zmetai = 6.295 have been clearly detected in the 
spectrum of GRB 050904. Therefore theoretical expecta¬ 
tion and possible ranges of zdla and zigm,u with respect 
to Zmetai should be discussed. 

The most clearly detected feature is SlI A1259.5 -|- Sill 
A1260.4, and Sll A1253.8 and Sill* A1264.7 are also de¬ 
tected beside the 1260 A feature. The profile fitting and 
curve of growth analysis revealed that these lines are not 
heavily saturated and only marginally resolved, giving 
column densities of logiVsn = 15.60^0 17 and logA^siii = 
14.29^039 and a velocity dispersion of ^ a few hundreds 
km/s (Paper I). The fine-structure silicon, Sill*, is often 
detected in GRB afterglow spectra (Vreeswijk et al. 2004; 
Ghen et al. 2005; Berger et al. 2005a, b) and indicates a 
high density environment, while it has never been clearly 
detected from DLAs in quasar (QSO) spectra. This sug¬ 
gests that the metal absorption system is a high density 
gas at the vicinity of the GRB explosion site (such as star 
forming regions, molecular clouds, or shells ejected by the 
GRB progenitor’s activity), rather than by absorption on 
a galactic scale of the host galaxy or by an intervening 
system. 

4 . 1 . Zdla versus Zmetai 

A clearly reasonable assumption is zdla = -^metai, be¬ 
cause no metal absorption lines were identified at red- 
shifts different from 6.295, except for the intervening sys¬ 
tem with a very different redshift, z = 4.840. Later we will 
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find that logiVni ~ 21.6 is necessary^ to explain the ob¬ 
served damping wing shape by the DLA absorption with 
■^DLA = -^metai- This column density is typical for DLAs 
found in GRB afterglows (Vreeswijk et al. 2004), and is 
also similar to those of molecular clouds in the Galaxy 
(Larson 1981; Solomon et al. 1987). This is consistent 
with the high density inferred from the Sill* detection. 
Adopting this hydrogen column density, we get the metal- 
licity estimates as [S/H] = —1.3 and [Si/H] = —2.9, using 
the solar abundance of Grevesse & Sauval (1998). Here 
we made no correction for ionization, assuming that the 
column densities of the singly ionized ions are equal to the 
total column densities, which is a reasonable assumption 
for silicon and sulphur in regions of large HI absorption 
like DLAs (Viegas 1995; Vladilo et al. 2001). Sulphur is 
not depleted to dust grains and the metallicity inferred 
from [S/H] is reasonable compared with those found in 
other GRB-DLAs. (See Vreeswijk et al. 2004 for a com¬ 
pilation of metallicities of GRB-DLAs and QSO-DLAs.) 
The large depletion of silicon compared with sulphur is 
not surprising, since large dust depletions compared with 
QSO-DLAs have been observed in a few GRB afterglows 
(Savaglio, Fall, & Fiore 2003). Therefore the observed col¬ 
umn densities and metallicities are reasonable in the case 
of the physical association between the DLA and absorp¬ 
tion lines. 

On the other hand, it seems rather unlikely that the 
DLA and the metal lines have different redshifts, again 
from the metallicity argument. In this case, no detectable 
lines associated with the DLA set 2a upper limits on 
metallicities as [S/H] < —1.5 and [Si/H] < —3.1, while the 
metallicity of the region associated with the metal lines 
would become larger than those estimated above. Hence 
there must be a large metallicity difference between the 
DLA and the metal absorption system. Such metallic¬ 
ity inhomogeniety seems rather unlikely if the DLA is 
also physically associated with the high density region 
in the vicinity of the GRB. On the other hand, if the 
DLA is a system on the galactic scale like QSO-DLAs 
in the host galaxy or in any intervening galaxies, we do 
not expect heavy depletion of Si (Savaglio et al. 2003), 
and hence the upper limit on [Si/H] should not be far 
from the total metallicity. However, such a low metallic¬ 
ity is not observed in any QSO-DLAs at z ~ 1-5, though 
there is no measurements of QSO-DLA metallicity at z > 6 
(Prochaska et al. 2003). Based on these considerations, 
we set zdla = -^metai in our baseline model. 

However, we conservatively allow the case of Zdla = 
■Zhost ^ -Zmetai, Supposing that the absorbing neutral hydro¬ 
gen is distributed on the scale of the host galaxy. Though 
a naive expectation is Zhost = ^metai, we expect a systemic 
velocity difference between them up to a few hundreds 
km/s depending on the host galaxy’s circular velocity 

^ This column density is slightly different from logA^jjl = 21.3 
adopted in Paper I, which is the best-fit DLA model when zdla 
is treated as a free parameter (see §5.2). The following metallic¬ 
ity estimates are correspondingly different from those in Paper 
I, but the difference of this magnitude does not affect the argu¬ 
ments in this paper. 


(l-^host — -^metaij $ 0,005). Furthermore, absorption lines 
blueshifted up to ~ 3,000 km/s with respect to the rest- 
frame of the host galaxy were observed in the afterglow 
spectrum of GRB 021004, which may have been acceler¬ 
ated by activities of the GRB or its progenitor (Mpller 
et al. 2002; Schaefer et al. 2003; Mirabal et al. 2003). 
Similar features were found also for GRB 020813 (Barth 
et al. 2003) and GRB 030226 (Klose et al. 2004). The 
blueshifted metal lines in the GRB 021004 spectrum were 
highly ionized (GiV and Siiv), but Lya absorption was 
observed as well at the same relative velocity. This indi¬ 
cates that the accelerated absorber is a mixture of neu¬ 
tral and highly ionized gas, and it is not unreasonable if 
blueshifted Si II is detected in some other GRBs. 

There are, however, two arguments disfavoring the 
blueshift of ~ 3,000 km/s. First, it seems rather unlikely 
that only metals in such a high velocity shell were de¬ 
tected without any detection at other velocities or the 
burst restframe. Second, such high velocity shells are ex¬ 
pected only within ~ 1 pc from the burster (Schaefer et al. 
2003; Mirabal et al. 2003), which might be inconsistent 
with the large amount of the silicon depletion, because 
strong GRB radiation would destruct dust grains within 
~ 1 pc of a burster (Waxman & Draine 2000; Fruchter 
et al. 2001; Draine & Hao 2002). However, a detailed 
modeling is required to examine these possibilities more 
quantitatively, and here we allow a blueshift of Zmetai up 
to (zhost — -^metai) ~ -fO.OT. Therefore, the possible range 
of Zdla becomes -0.005 < (zdla - ^^metai) $ +0.07. 

4-2. ZIGM.U versus Zmetal 

We should set zigm.u = -Zhost if the GRB occurred in a 
galaxy embedded in the IGM that had not been ionized 
by the emission from the galaxy. Barkana & Loeb (2004) 
estimated that typical z ~ 7 GRBs occur in galaxies whose 
dark halo mass is Mh ~ 4 x IO^Mq, based on the standard 
structure formation theory. They found that the proper 
(i.e., not comoving) radius of the ionized bubble created 
by stellar radiation from such galaxies would be $ 0.1 
Mpc, corresponding to a wavelength shift of AAobs ~ 3A or 
(•ziGM.u — -Zhost) ~ —0.0025. Instead, infall of IGM into the 
host galaxy may result in a relative redshift of the IGM ab¬ 
sorption, with a wavelength shift of AAobs $ lA for a typi¬ 
cal GRB host galaxy (Barkana & Loeb 2004). Since these 
are sufficiently small compared with the wavelength range 
of the damping wing, we set Zigm.u = Zhost in the base¬ 
line model. On the other hand, Haiman (2002) estimated 
the proper size of ionized bubbles around LAEs detected 
in deep surveys as ^ 0.8 Mpc (zigm,u — -Zhost ~ —0.02), 
assuming a star formation rate (SFR) of ~ lOMQ/yr (as 
typically estimated for LAEs from the Lya emission lumi¬ 
nosity) and stellar age of 10® yr. This means a stellar mass 
greater than M* ~ 10®M@, which is considerably larger 
than that estimated by Barkana & Loeb (2004). Hence 
we consider that such a large ionized bubble is rather un¬ 
likely for typical GRBs at z ~ 7. It should be noted that 
the LAEs found in deep surveys are not necessarily typical 
galaxies or a tracer of typical star formation at z ~ 7. 

Based on these considerations combined with the pos- 
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sible difference of Zhost from Zmetai discussed in §4.1, we 
define the possible range of zigm.u as —0.02 ^ (zigm.u — 
-^metai) ^ +0.07, with a canonical value of zigm.u = -^metai- 

5. Fitting to the Red Damping Wing 

5.1. Fitting Procedures 

In this section we will fit the models described above to 
the observed shape of the red damping wing. The fit is 
performed in a wavelength range of Aobs = 8925-9938A. 
The lower bound corresponds to the wavelength where 
the wing drops sharply to the flux zero level. The wave¬ 
lengths shorter than this are excluded since the noise by 
atmospheric lines is relatively high and there may be Lya 
emission from the host galaxy (see discussion in §7.2). We 
then excluded the following wavelength regions of known 
absorption lines: 9033-9046A (likely to be CIV in an inter¬ 
vening system at z = 4.840), 9137-9158A (Sll A1253.8), 
9180~9206A (Sii A1259.5 + Sill A1260.4), 9218-9238A 
(Sill* A1264.7), 9490-9505A (Ol A1302.1), and 9728- 
9750A (Cii A1334.5). In addition, we excluded a range 
of 9305-9490A, since there are some absorption features 
that are probably atmospheric. There are 268 data points 
(pixels) in the final wavelength ranges used in the anal¬ 
yses, which are indicated by thick horizontal lines in Fig. 
1 . 

Then the model spectrum is fit to the data, with nor¬ 
malization always chosen to minimize the value. In 
addition to the DLA and IGM absorptions, the Galactic 
extinction by E{B — V) = Q.QQQ mag (Schlegel, Finkbeiner, 
& Davis 1998) is also taken into account in the model 
spectrum. The spectral resolution is sufficiently narrower 
than the observed damping wing shape, and hence we 
ignore the smoothing effect by it. Figure 3 shows the his¬ 
togram of the deviation of the data flux in a pixel from 
the model flux, (Fdata - 7Aiodei)/Aerror, for One of the ac¬ 
ceptable fits with zdla = 6.295, logAini = 21.62, and 
a:Hi = 0 (x^ = 277.29). The distribution can well be fit 
by the Gaussian distribution, indicating that the error is 
well controlled and we can derive a reasonable statistical 
constraint on the model parameters by the standard x^ 
analyses. 

5.2. Separate Fittings by DLA and IGM Absorption 

As the first simple analysis, we try to fit the two com¬ 
ponents separately to the data. Figure 4 shows the confi¬ 
dence region of iVHi--ZDLA foi' the DLA fit and a:Hi-2iGM.u 
for the IGM fit. The best fits are (logAini, z^dla) = 
(21.34, 6.320) with x^ = 265.54 and (xhi, -sigm.u) = 
(1.58, 6.348) with x^ = 270.79, respectively. 

The best-fit value of Zdla is different from Zmetai by 
~ +0.025, which is considerably larger than that ex¬ 
pected from the circular velocity of the host galaxy, in¬ 
dicating a possibility that the metal absorption system 
was accelerated by the GRB or its progenitor. However, 
Zdla = ^metai = 6.295 is also marginally consistent within 
the confidence region of 99.7% G.L., and taking into ac¬ 
count the arguments given in §4.1, we keep this case as 
the baseline model. 


Interestingly, almost neural IGM with xm ~ 1 extend¬ 
ing to ZiGM,u ~ 6.36 gives a marginally acceptable fit for 
the red damping wing, indicating a possibility that we 
may have observed for the first time the signature of the 
completely neutral IGM before the reionization. The red- 
shift difference of (zigm.u — zmetai) ~ +0.065 is within the 
theoretically possible range discussed in §4.2. 

We show the model spectra for the two representative 
cases: “the DLA-dominated model” with zdla = 6.295 
and (logA^m, ^igm.u, a;Hi) = (21.62, 6.295, 0) in Fig. 1 
and “the IGM-dominated model” with (—oo, 6.36, 1.0) 
in Fig. 2. Although the red damping wing shape can 
be reproduced even with xhi = 0 in the DLA-dominated 
model, the IGM absorption with xm ^ 10“^ and zjgm.u ^ 
6.26 are required for the absorption of photons blueward 
of the Lya resonance. 

5.3. DLA-IGM Joint Fitting 

Next we try a joint fit including both the DLA and IGM 
absorptions. Here we treat zigm.u as a free parameter but 
assume zdla = ^metai = 6.295 to reduce the number of 
free papameters. Figure 5 shows the allowed region for 
A^hi and xhi, for several values of zigm.u- The best-fit is 
obtained at (zigm.u, log-Nni, xm) = (6.36, 21.59, 0.016) 
with x^ = 269.70. Gonstraining xhi > 10“^ as inferred 
from the GP trough in quasar spectra at z ^ 6, we obtain 
an upper limit of zigm.u + 6.371 at 99% G.L., otherwise 
the IGM absorption will erase the light transmission at 
Aobs > 8950 A. 

As expected from the single component fits, the damp¬ 
ing wing can be explained mostly by the DLA absorption 
at Zdla = 6.295 with negligible IGM absorption, while it 
can also be fit dominantly by the IGM absorption with 
xhi ~ 1 and Zigm.u ~ 6.36. These results tell us that there 
is an unfavorable degeneracy between absorptions by a 
DLA and IGM (Miralda-Escude 1998; Barkana & Loeb 
2004). However, in the next section we show that this 
degeneracy can be broken by utilizing the features around 
the Ly/3 resonance. 

5.4-. Constraints from the Ly(3 Feature 

Here we apply the same absorption models of the DLA 
and IGM, but for the Ly/3 transition at the restframe 
wavelength of Xp = 1025.72 A. We use the same formulae 
given in §3, but replace A^ and fa by the corresponding 
values for Ly/3. The Ly/3 absorption oscillator strength is 
fp = 0.07910 and gu/gi = 3 is the same for all the Lyman 
series transitions by the selection rules (Gox 2000; Morton 
1991). The frequency-integrated cross section f cr(z/)di/ is 
reduced by fp/fa = 0.190, and the width of the damping 
wing in the wavelength space (cx ApX^) is reduced by the 
same factor. The optical depth of the red damping wing 
("Tigm) is reduced by a factor of {fpjfaff = 0.0361 for the 
same Zobs (= Aobs/A/3 - 1 for Ly/3). 

The predictions for the Ly/3 absorption feature are 
presented correspondingly to Lya in Fig. 2. The fea¬ 
ture of the DLA-dominated model is in good agreement 
with the data, while the emission (or, residual of absorp¬ 
tion) feature around 7500-7520A is clearly inconsistent 
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with the IGM-dominated model, in which the transmis¬ 
sion is completely attenuated at wavelengths blueward of 
A/3(1 -I- ziGM.u) ~ 7550A. By visual inspection of the orig¬ 
inal CCD image of the spectrum (shown in Fig. 2), we 
confirmed that the feature of 7500-7520A is the real sig¬ 
nal. Therefore the Ly/3 feature safely excludes the IGM- 
dominated model. By requiring that the IGM Ly/3 ab¬ 
sorption does not erase the feature at 7500-7520A, we set 
a firm upper limit of zigm.u A 6.314, and within this con¬ 
straint, the IGM absorption has only small contribution 
to the observed red damping wing even if a;Hi = 1 (see 
dashed curves of Fig. 2). However, this contribution is 
not negligible, and in fact we can derive some constraint 
on xhIj as will be shown below. 

6. Implications for the IGM Neutral Fraction 
around GRB 050904 

After the long analyses presented above, we finally ar¬ 
rived at the most likely interpretation of the red damping 
wing: it is mainly contributed by the DLA associated with 
the metal absorption lines at Zmetai = 6.295, and the contri¬ 
bution from neutral IGM is not necessary at all. However, 
as shown in Fig. 2, if the IGM is completely neutral, 
the IGM absorption with zigm,u ~ 6.295 has small but 
non-negligible contribution to the observed damping wing 
shape. Then the next question is which value of xhi is fa¬ 
vored, under the constraints of zigm,u = -^dla = 6.295 but 
A^hi is treated as an unknown free parameter and always 
chosen to minimize for each value of xhi- 

For parameters other than xm and A^hi, we take 
^^iGM.u = ^DLA = 6.295, ziGM.i = 6, Po = “1.25, and Ay = 0 
as the baseline model parameters. We then find that 
becomes minimum at xm = 0 with logA^ni = 21.62, and 
it continuously increases with xhi, changing the best-fit 
value of logA^Hi to 21.54 at xhi = 1- [The preference to 
smaller xm for zigm.u ^ 6.32 was also found in the results 
of the DLA-IGM joint fit (see §5.3 and Fig. 5).] Then 
formally we can estimate the confidence region of xm by 
requiring that /:^x^{xm) = X^(a;Hi) - xiin obeys the chi- 
square distribution with one degree of freedom (Lampton, 
Margon, & Bowyer 1976; Press et al. 1992). Table 1 
presents the x? values, the exclusion confidence level of 
the xm = 1 model corresponding to Ax^(l), and upper 
limits on xhi with 68 (Icr), 95, and 99 % C.L. We find the 
best fit value of a;Hi = 0.00 in all models and upper limits 
of Xhi <0.17 and 0.60 (68 and 95% G.L., respectively) for 
the baseline model. 

In order to show why the a;Hi = 1 model gives a worse 
fit, we plot in Fig. 6 the model spectra and the observed 
data in the form of the residual from the best-fit model 
of xhi = 0 (the DLA-dominated model). The worse fit of 
the xm = 1 model comes from the combination of the two 
effects: the change of the damping wing shape in Aobs = 
8900-9200 A is not favored from the data, and the model 
predicts higher flux than the xhi = 0 model at the longest 
wavelength range of Aobs ^ 9500A, while the data favor 
even lower flux than that of the xhi = 0 model. 

Before we conclude that the case of a;Hi = 1 is not fa¬ 


vored, however, any systematic uncertainties must be ex¬ 
plored. The results of this uncertainty check will be given 
below. Table 1 shows the fit results changing one or two 
of the model parameters while the other parameters are 
kept to the values of the baseline model. (The DLA col¬ 
umn density Nm is treated always as a free parameter 
minimizing y^.) We found that the x^ minimized at 
xm = 0 in all cases, favoring smaller values. 

6.1. Spectral Index of the Afterglow SED 

The change of the afterglow spectral index, /3o, within 
the quoted error of A/3obs = ±0.25 (Tagliaferri et al. 2005) 
does not have a significant effect. As can be seen in Fig. 6, 
the data favors a bluer spectrum than the baseline model 
with /3o = —1.25, and hence a larger /3o gives a better fit. 
However, this effect is not large; the change of A/3o = 0.25 
in the wavelength range of 9000-10000A results in the 
fractional change of ~ 2.4% in F\, or ~ 0.024 in Fig. 6. 
It should be noted that the quoted error of Haislip et 
al. (2005) is even smaller (A/3obs =l'lo i 4 )- If the cooling 
frequency break in the SED passed through the NIR bands 
at an epoch between the photometric observations and 
our spectroscopic observation, the afterglow SED would 
become even redder, which is not favored from our data. 

6.2. Dust Extinction in the Host Galaxy 

Though the observed NIR colors do not show evidence 
for extinction by dust in the host galaxy (Tagliaferri et al. 
2005), the wavelength coverage by photometric observa¬ 
tions redward of the Lyman break is rather narrow, and 
substantial extinction may still be possible. Therefore we 
tried a model with substantial extinction, where the orig¬ 
inal afterglow spectrum jSo is chosen so that the expected 
NIR colors are consistent with the observed ones. The 
change of observed {J — K) color at z ~ 6.3 can be related 
to the restframe V magnitude extinction as A( J — K) = 
0.68Ay and 1.83Ay for the extinction curves of the Milky 
Way (MW) and the SMG, respectively (e.g., Pei 1992). 
On the other hand, the change of the afterglow spectral 
index is related to the NIR color as A{J — K) = — 0.61A/3o. 
Therefore the expected (J — K) color is not changed if 
we introduce extinction with changing /3o by a relation. 
Ay = 0.90A/3o (MW) or 0.33A/3o (SMC). In Table 1, we 
show the fit results by models with Ay = 0.45 (MW) and 
0.17 (SMG) corresponding to A/3o = ±0.5 (/3o = —0.75). 
Here we have used the functional form of extinction A(A) 
around A ~ Aq given in Barkana & Loeb (2004). We found 
that both x^ and Ay^ are increased by introducing extinc¬ 
tion by dust, indicating that substantial extinction is not 
favored and does not change the preference to lower xhi 
values. 

6.3. The Redshift Parameters 

Ghanging zigm.u within the possible range of 6.27- 
6.314 does not change the preference to lower xbi values. 
Because Zigm,i is rather uncertain, changing this param¬ 
eter has a relatively large effect. Larger values of zigm,i 
result in a smaller amount of the neutral IGM and hence 
the upper limit on xhi becomes weaker. However, the 
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completely neutral IGM can still be excluded by ^^90% 
C.L. even for zigm,i = 6.2. 

To show the characteristic redshift range of IGM neutral 
hydrogen contributing to tigm , we plot tjgm as a function 
of ^iGM.i with ziGM,u = 6.295, for three values of Aobs = 
8925 (the lower bound of the wavelength range used in the 
fitting analysis), 8970, and 9050 A, in Fig. 7. The IGM 
optical depth is shown as a ratio to that with zigm.i = 0, 
which is Tigm = 0.37, 0.20, and 0.10 for the three values 
of Aobs, respectively. From this figure we see that more 
than half and 80% of the IGM optical depth is contributed 
from hydrogens at z ^ 6.2 and 6.0, respectively. 

As a check for the possibility that zdla is different 
from Zmetai but Zdla = -^host, we tested the models with 
zdla = -ziGM,u = 6.29 and 6.314. This range can be de¬ 
rived assuming Zhost = ^igm.u and taking into account the 
possible ranges of Zhost and zjgm.u with respect to Zmetai, 
as well as the constraint from the transmission around the 
Ly/3 wavelength. We found that the preference to low xhi 
is not changed. 

6.4- Weak Unidentified Absorption Lines 

It is difficult to quantitatively estimate the possible sys¬ 
tematic effects by weak absorption lines that cannot be 
discriminated from noise. However, if such weak lines are 
distributed uniformly in the htting wavelength range, the 
effect is simply changing the overall normalization that 
would not affect our conclusions. We also repeated the 
same calculation without removing the wavelength ranges 
of the identified absorption feature, to see the sensitivity 
to the absorption features (see the “Lines Included” row 
of Table 1). The data points are now 379, and is un¬ 
acceptably large due to the absorption features that are 
not included in the model. However, we find a similar 
Ax^(l), and hence the preference to small xhi values is 
not sensitive to the treatment of absorption lines. 

The absorption by vibrationally excited molecular hy¬ 
drogen may affect the spectrum around the restframe Lya 
(Draine 2000; Draine & Hao 2002). In fact, Haislip et 
al. (2006)^ proposed that the anomalously attenuated 
Z-band flux by a factor of ~ 3 of the early afterglow of 
GRB 050904 is due to the molecular hydrogen absorp¬ 
tion. However, this signature completely disappeared by 
~ 3 days, and hence it is unlikely to affect our spectrum 
taken 3.4 days after the burst. The molecular hydrogen 
signatures have not yet been clearly detected in other af¬ 
terglow spectra (e.g., Schaefer et al. 2003; Vreeswijk et 
al. 2004). 

6 .5. Time Variability 

Ionization by strong afterglow flux may result in the 
time variability of the amount of hydrogen and metal 
absorptions (Perna & Loeb 1998; Draine & Hao 2002). 
Though statistically significant variability was not found 
in the observed spectrum of GRB 050904, variability 
within the errors cannot be excluded. If such a variabil¬ 
ity of Ahi is present, the observed shape of the damping 
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wing is a superposition of the shapes with different val¬ 
ues of Ahi- It may induce systematic bias in the analysis 
presented here, since the damping wing shape depends 
non-linearly on iVni. As a test of this possibility, we tried 
an extreme model where the spectrum is a mean of the 
two spectra with different values of the column density, 
N+ and N -, as log N± = log Ari ± A log Ari • We take 
AlogARi = 0.2 as the variability allowed within the ob¬ 
servational error, which is inferred from the statistical er¬ 
ror of logARi obtained by the fitting analysis (see Fig. 
4). We then find that the preference to lower xri is 
hardly changed even for this rather extreme model (see 
the “Variability Check” row of Table 1). 

It should be noted that the neutral hydrogen shell must 
be very close to the burster, as r ^ 4 x 10^^ cm, in order for 
the column density of Ari ~ 21.6 to substantially change 
during the spectroscopic observation by the ionizing flux 
estimated from the observed Fx. This scale is close to the 
location of the external shock at 3.4 days in observer’s 
time for typical afterglow parameters and z = 6.3 (e.g.. 
Sari, Piran, & Narayan 1998). We derived in Paper I an 
estimate of electron density rie ~ cm“^ from the 

observed hne-structure Si* line, assuming collisional exci¬ 
tation by electrons. The ionization fraction is unknown, 
and if ionization fraction is much smaller as ne/nn $ 10“^, 
the excitation by neutral hydrogen becomes more impor¬ 
tant, leading to a density estimate of nn 10^ cm ^ (Silva 
& Viegas 2002). Combined with the DLA column den¬ 
sity, the DLA scale of ~ 0.1-1 pc is inferred. On the 
other hand, Berger et al. (2005b) suggested a possibility 
of fine-structure excitation by GRB afterglow radiation 
for the case of GRB 051111, deriving a larger distance 
scale of ~10-20 pc. Therefore the scale estimate is rather 
uncertain, but it seems that there is a parameter space 
of the DLA where no Ari variability is expected but it is 
close enough to be consistent with the Sill* detection. 

Based on these results, we conclude that there is an 
evidence for considerable ionization of the IGM around 
the host galaxy of GRB 050904 already at z = 6.3. 

7. Discussion 

7.1. Prediction for the Case of Low Ari 

We found that the DLA column density is log Ari > 21 
for GRB 050904, and in such a case the DLA absorption 
dominates the absorption by neutral IGM extending to 
the same redshift as that of the DLA. Therefore it is dif- 
hcult to derive a very strong constraint on a;Hi, though 
we obtained some constraints as presented above by de¬ 
tailed statistical analyses. Vreeswijk et al. (2004) com¬ 
piled seven GRBs having estimates of Ari from afterglow 
spectra, and though the majority of them have column 
densities larger than log Ari > 21, two (GRB 011211 and 
GRB 021004) have low column densities of log Ari ^ 20. 
If we detect such a GRB at z ^ 6, it would provide an op¬ 
portunity to probe the reionization with a much cleaner 
environment. 

To demonstrate this possibility, we show a model with 
(log ArrXri) = (20.0,1.0) with ZdlA = ^IGM.u = ^^metal = 
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6.295, in Fig. 8. In this case the shape of the damp¬ 
ing wing is dominantly contributed by the neutral IGM 
absorption. Also shown is DLA absorption at the same 
redshift but with various values of A^hi- It can be seen 
that the shape of IGM absorption is different from that 
of DLAs with any values of A^hIi making it possible to 
prove the dominance of the IGM component. Another 
possibility to explain such a data by DLA absorption is 
to consider a blueshifted absorption by a DLA with a rel¬ 
ative velocity of a few thousands km/s with respect to 
-^metai, as we observed the degeneracy between the DLA 
and IGM absorptions in §5.2. Such a velocity may again 
be achieved by acceleration of absorbing clouds by the ac¬ 
tivity of GRBs or their progenitors, but such absorbers are 
likely to be polluted by metals, and we expect blueshifted 
metal lines as well, in addition to those in the restframe 
of the host galaxy. Therefore we can discriminate this 
case as well, and hence high-quality spectra of 1ow-7Vhi 
afterglows in the future will provide us with an invaluable 
opportunity to measure xm more accurately. 

1.2. Lya Emission from Host Galaxies and Reionization 

We expect Lya emission in the spectrum, if the huge 
Hi absorption is by gas in the vicinity of the GRB and 
the host galaxy has detectable Lya emission like GRB 
030323 (Vreeswijk et al. 2004). The spectrum of GRB 
050904 does not show a clear Lya emission (§2) and we 
set an upper limit on the Lya flux of F(Lya) < 1.7 x 
10“^® erg cm“^s“^ for a width of 8.5 A, which is the 
FWHM of the spectral resolution and corresponding to a 
velocity dispersion of ^ 300 km/s as well. Then the upper 
bound on the Lya luminosity becomes L(Lya) < 7.9 x 10^^ 
erg/s at z = 6.295, which can be translated into extinction- 
uncorrected SFR of < O.79M0/yr, by using the empirical 
relation (Kennicutt 1998; Cowie & Hu 1998). 

It has been indicated that the Lya emission from GRB 
host galaxies is stronger in terms of the equivalent width 
(EW) than that from Lyman-break galaxies (LBGs) at 
similar redshift (z 3); all GRB host galaxies at z ^ 2 
are consistent with having EWrest(Lya) > 10 A (Fynbo 
et al. 2003; Jakobsson et al. 2005), which is as large as 
that of LAEs found in deep surveys. For comparison, only 
about a third of LBGs at similar redshifts have such large 
EW (Shapley et al. 2003). Then it is suggested that Lya 
emission from host galaxies of GRBs at z ^ 6 may also 
give some information on the reionization. 

To demonstrate this possibility, we show the IGM ab¬ 
sorption profile with zjqm.u = 6.295 for xm = 0.01,0.2, and 
1.0 in Fig. 2. If xhi ~ 1 and the Lya line center is the same 
as Aq,(1-|-zigm,u), the Lya photons within ^10 A from the 
line center are heavily absorbed. For comparison, LAEs 
detected at z ~ 6 show line widths of ~ lOA, corresponding 
to a few hundreds km/s (Kodaira et al. 2003; Taniguchi 
et al. 2005). Even larger velocity dispersion may also 
be possible; a redward tail of Lya emission extending be¬ 
yond a relative velocity difference of ^ 1000 km/s has been 
found in a starburst galaxy at z ~ 3 (Pettini et al. 2000). 
Then we expect that, if Ihi 1, the Lya line luminosity 
is strongly attenuated allowing transmission only in the 


redward tail. On the other hand, if xm ^ O.OI, about half 
of the Lya emission is transmitted. Therefore statistical 
comparison of Lya line emissivity of GRB host galaxies 
at z ^ 6 and ^ 6 may give an interesting information for 
the reionization. 

As a reionization probe based on Lya emission, GRB 
host galaxies have an advantage of no selection effect com¬ 
pared with LAEs selected by Lya emission. Brightness of 
afterglows allows us to accurately measure the redshift by 
metal absorption lines, as demonstrated by GRB 050904, 
and hence the plausible Lya line center can be deter¬ 
mined more accurately than LAEs or LBGs. It should 
be noted that, even though LAEs at z 6 have asym¬ 
metric Lya emission profiles with a sharp cut-off at the 
blueward side, a strong conclusion about the IGM neutral 
fraction cannot be derived because of the general lack of 
information about the line center (Haiman 2002). On the 
other hand, since the brightness of host galaxies is not 
relevant to the detectability of GRBs, typical GRBs may 
occur in smaller galaxies than LAEs/LBGs, and hence 
typical absolute Lya luminosity may be smaller, in spite 
of the suggestion of large EW for GRB host galaxies. (See 
§4.2 for a discussion about the typical mass of GRB host 
galaxies at z ^ 6 expected from the structure formation 
theory.) A larger sample of GRB host galaxies at z ^ 6 is 
necessary to investigate more about these possibilities. 

Concerning the case of GRB 050904, a search for 
the host galaxy in the z' band will be interesting 
in this context. Assuming EWrest(Lya) > 10 A as 
suggested for GRB host galaxies, the obtained upper 
bound on Lya emission requires that the continuum 
level of the host galaxy must be lower than F\ < 2.3 x 
I0“^°(EWrest/10A)“^ erg cm“^s“^A“^. Since the Lyman 
break is in the midst of the z' band filter, we take into ac¬ 
count a dimming by a factor of about 2, and we get the 
corresponding AB magnitude of z' > 27.7. Therefore, if 
follow-up observations find the host galaxy of GRB 050904 
brighter than z' = 27.7, it would mean that EWrest(Lya) < 
lOA, which is smaller than those found in GRB host galax¬ 
ies at z ~ 3. Thus it indicates a possibility of attenuation 
by neutral IGM, though only one case is not sufficient 
and statistical studies are required. Galaxies brighter than 
z' = 27.7 can be detected by imaging observations of exist¬ 
ing large telescopes; typical magnitude of LAEs at z 6.6 
found in the Subaru Deep Field is z' ~ 26-27.8 (Taniguchi 
et al. 2005). 

8. Summary and Conclusions 

We presented a comprehensive theoretical modeling of 
the red damping wing of the Lya absorption found in 
the optical afterglow spectrum of GRB 050904 at z 6.3, 
which provides the first opportunity of studying the cos¬ 
mic reionization by using GRBs. We tried to model the 
observed damping wing shape by the two components of 
absorbers: one is by the DLA associated to the GRB host 
galaxy, and the other is by neutral hydrogen in the IGM. 
The redshift of the metal absorption lines in the spectrum 
is z^metai = 6.295 ± 0.002, but we allowed different values of 
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zdla (the DLA redshift) and 2 igm,u (the upper extension 
bound of neutral hydrogens in the IGM), and discussed 
various theoretical possibilities for the deviation of these 
parameters from Zmetai- 

The shape of the red damping wing can be explained 
either by the DLA at zdla ~ -^metai with logA^ni ~ 21.6, 
or by almost neutral IGM extending to a higher redshift 
of ziGM.u ~ 6.36. Though the DLA seems a more straight¬ 
forward interpretation, we cannot exclude the latter pos¬ 
sibility simply by the redshift difference, since blueshift of 
metal absorption lines up to a few thousands km/s with 
respect to the restframe of the host galaxy has been ob¬ 
served in a few GRBs. 

However, we found that the Ly/3 feature can be used 
to break this degeneracy, since the two different solutions 
predict different wavelengths at which the Ly/3 GP trough 
ends. Then we concluded that the damping wing is mostly 
contributed from the DLA at zdla ~ -^metai, and derived a 
firm upper bound of zigm.u < 6.314. We argued that the 
DLA is likely to be associated physically with the metal 
absorption lines, and the inferred column densities, metal- 
licities, and depletion of silicon are all reasonable as a DLA 
found in a GRB afterglow. 

Next we examined the preferred value of the IGM neu¬ 
tral fraction, xm, in the viable model of zigm.u = ^dla = 
6.295. Treating A^hi of the DLA as a free parameter, we 
found that a smaller value of cchi is favored with the best- 
fit value of xhi = 0.00, and upper limits of cchi < 0.17 and 
0.60 (68 and 95% G.L., respectively) were derived. We 
examined various possible systematic uncertainties that 
could affect this result, including the afterglow spectral 
index, dust extinction at the host galaxy, the redshift pa¬ 
rameters of the DLA and IGM absorptions, weak uniden¬ 
tified absorption lines, and time variability of the DLA col¬ 
umn density. We found that none of these effects is large 
enough to change the above result. Hence we conclude 
that the universe was largely ionized already at z ~ 6.3, 
excluding the completely neutral IGM at ~ 99% C.L. 

This is the first quantitative upper limit on xhi at z ^ 6 
by a direct method,^ being consistent with the recent re¬ 
sults by an indirect approach using the number density 
evolution of LAEs (Malhotra & Rhoads 2004; Stern et 
al. 2005; Haiman & Gen 2005). Since the IGM is op¬ 
tically thin for photons in the damping wing region, all 
the IGM neutral hydrogens at z ~ 6.1-6.3 (see Fig. 7) 
contribute to the damping wing, allowing us to derive a 
robust constraint on the mass-weighted xhi which is in¬ 
sensitive to any dumpiness of IGM within this redshift 
interval. Combined with the suggestions of cchi ^0.1 from 
quasar spectra (Mesinger & Haiman 2004; Wyithe, Loeb, 
& Carilli 2005), a plausible value of xhi ~ 0.1 is suggested 
for the IGM at z ~ 6-6.3. 

Though there are some “data points” of the IGM optical depth 
at 2 ^ 6 derived directly from the GP trough of quasar spec¬ 
tra (Fan et al. 2006), these are by averaging sharp spikes of 
transmission in a wavelength range, which are probably corre¬ 
sponding to small regions of ionized bubbles. Therefore there is 
no upper limit on the mass-weighted or volume-weighted optical 
depth. 
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The large DLA column density of logA^ni ^21 domi¬ 
nates the IGM absorption extending to the same redshift, 
making it difficult to derive a stronger constraint on cchi 
than derived here. However, some GRBs have low column 
densities of logiVni ^ 20, and detection of such GRBs at 
z ^ 6 will be a promising chance to get better information 
for the reionization history of the universe. 

We did not detect Lya emission from the host galaxy, 
leading to an upper limit for extinction-uncorrected star 
formation rate as SFR < O.79M0/yr. We discussed the 
potential of Lya emission from GRB host galaxies as a 
reionization probe. Statistically smaller equivalent width 
and transmission only at the red tail of Lya emission are 
expected for GRB host galaxies before the reionization, 
compared with those at lower redshifts. This may be 
tested by using future large samples of GRBs at z > 6. 
As a reionization probe using Lya emission, GRB host 
galaxies have an advantage of being free from the selection 
bias compared with LAEs. Another advantage compared 
with both LAEs and LBGs is that an accurate redshift 
determination is possible by absorption lines in bright af¬ 
terglows. On the other hand, a possible disadvantage is 
that typical GRB host galaxies and their absolute Lya lu¬ 
minosity may not be as bright as LAEs and LBGs found 
in deep surveys. 

We conclude that the GRB 050904 has opened a new era 
of cosmological study by GRBs, and future data will give 
us even more unique and important information about 
the epoch when the early-generation luminous objects 
changed the physical state of almost all the baryonic mat¬ 
ter in the cosmos. 
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Table 1. Constraint on xhi for Various Models 


Models 

Amin 

X^(a;Hi = 1) 

= 1) 

C.L.(%)t 

upper limits on xm* 

68 95 99 C.L. (%) 

Baseline^ 

277.29 

284.12 

6.84 

99.1 

0.17 

0.60 

0.98 

/3o = -l/-1.5 

273.86/281.10 

280.07/288.75 

6.21/7.65 

98.7/99.4 

0.18/0.15 

0.66/0.52 

1.08/0.88 

Av = 0.45/0.17^ 

316.93/283.18 

329.72/291.12 

12.80/7.94 

99.97/99.5 

0.09/0.15 

0.32/0.52 

0.54/0.85 

aiGM.u = 6.27/6.314 

277.29/277.29 

283.55/283.21 

6.27/5.93 

98.8/98.5 

0.19/0.19 

0.63/0.68 

1.06/1.12 

2iGM,i = 5.5/6.2 

277.29/277.29 

286.70/280.07 

9.41/2.78 

99.8/90.5 

0.13/0.38 

0.45/1.36 

0.74/2.22 

ZDLA = 6.29/6.314" 

281.47/266.41 

288.49/271.25 

7.03/4.83 

99.2/97.2 

0.16/0.27 

0.57/0.84 

0.96/1.27 

Lines Included 

796.03 

802.53 

6.50 

98.9 

0.17 

0.63 

1.02 

Variability Check 

283.98 

289.54 

5.56 

98.2 

0.20 

0.71 

1.18 


*See §3 for the reason why the apparently unphysical cases of x^n > 1 are allowed here. 

^In all models, increases with x^n and the minimum is realized at xjji = 0. 

^The exclusion confidence level for the case of xhi = 1- 

^The baseline model parameters are: 2 :igm,u = -s^DLA = 6.295, ^^igm 1 = 6, /5o = —1.25, and Ay = 0. The second row and below show the 
models when some parameters are changed from the baseline model. 

^The spectral index is changed into /3o = —0.75 to keep the expected NIR colors consistent with the observed ones for the MW/SMC 
extinction curves, respectively. 

II The IGM redshift parameter 2 :igm,u is kept to be the same with 2 :dla- 



r T— T- ■ -i^~i-r —■ I ■ I ■ I ■ I 

6600 6600 9000 9200 9400 0600 9600 lOOOO 


Observed Wavelength [A] 


Fig. 1. The afterglow spectrum of GRB 050904 taken 3.4 days after the burst. The spectrum is binned by 3 pixels, and la errors 
are also shown with an offset of —1.0. The GCD image of the spectrum in the corresponding wavelength range is shown at the 
top of the figure. The Lya resonance and identified absorption lines are indicated with the redshift ^metal = 6.295, except for the 
intervening CiV system at 2 = 4.840. The thick horizontal lines in the upper right region show the wavelength ranges used in the 
spectral fitting, where identified absorption features are removed. The solid curve shows the model absorption by a DLA with 
logNm = 21.62 and 2 :dla = 6.295 (the DLA-dominated model). The dotted line shows the original unabsorbed spectrum of the 
afterglow, with the spectral index of 00 = —1.25. The dashed curve shows the model absorption by the IGM with 2 igm u = 6.295 
and xhi = 10~^, which is almost a vertical line at the Lya resonance. The Galactic extinction of E(B — V) = 0.060 mag is taken 
into account in all the model curves. 
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Fig. 2. The same as Fig. 1, but the close-up of the Lya (upper panel) and Ly/3 (lower panel) regions. The spectrum is binned 
by 2 pixels, and Icr errors are shown by dashed lines with an offset of —0.5. The wavelength scales of the two panels are chosen so 
that the Lyo: and Ly/3 wavelengths at a fixed redshift coincide with each other on the horizontal scale (i.e., Atop/Aiow = Aa/A^). 
The vertical solid line marks the Lya and Ly/3 resonances at z = 6.295. The CCD images of the spectrum in the corresponding 
wavelength ranges are shown at the top and bottom of the figure. The solid and dotted curves are the same with those in Fig. 1. 
The three dashed curves show the model absorption by the IGM (2 igm,u = 6.295), with = 0.01, 0.2, and 1.0, as indicated. The 
dot-dashed curve is the model absorption by the IGM but with 2 :igm,u = 6.36 and a:Hi = 1-0 (the IGM-dominated model). 



(F 


data 


:.) / F 


error 


Fig. 3. The histogram for the deviation of the observed flux in a pixel from a model spectrum, normalized by the Icr error of the 
data. The reference model is with logA^Hi = 21.62, 2 :dla = 6.295, and thi = 0 (the DLA-dominated model). There are 268 data 
points (corresponding to the number of pixels) in the wavelength range used in the fitting analysis. The solid curve is the fit by the 
Gaussian distribution with cr = 1. 
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Fig. 4. The arrowed regions for the separate fits of the DLA (upper panel) and IGM (lower panel) absorption models to the red 
damping wing. The horizontal axis is the DLA redshift (^dla) foi" DLA model but is the upper redshift bound of the extension 
of neutral hydrogens in the IGM (2 :igm,u) for the IGM model. The contours show the levels of = 2.30,9.21, and 18.4 from the 
minimum corresponding to 68.3 (dotted), 99 (thick solid), and 99.99% (thin solid) confidence levels for two degrees of freedom. 
See text (§3) for the reason why the apparently unphysical region of xhi > 1 is also shown. The redshift of metal absorption system 
and its error are indicated by the thick and thin vertical lines, respectively. 
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Fig. 5. The arrowed regions for the joint fit of the DLA and IGM absorption models. The four panels are for different values of 2 :igm,u 
as indicated in the panels. The DLA redshift is fixed to 2 :dla = -Smetai = 6.295. The contours show the levels of Ax^ = 3.53,11.3, 
and 21.1 from the minimum x^? corresponding to 95 (dotted), 99 (thick solid), and 99.99% (thin solid) confidence levels for three 
degrees of freedom. See text (§3) for the reason why the apparently unphysical region of xhi > 1 is shown. 
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Fig. 6. The model and observed spectra of GRB 050904, in the form of the residual from the DLA-dominated model 
[(^dlAj log■2^IGM u! ^Hi) = (6.295, 21.62, 6.295, 0)]. The connected pixels are binned by up to 20 pixels, after removing 
the wavelength ranges which were not used in the analyses because of absorption features. The thick solid curve is the best-fit 
model (logA^Hi = 21.54) under the constraints of xhi = 1 and zyu^a. = ^iGM,u = 6.295. 



Fig. 7. The IGM optical depth as a function of ^igm,!? which is normalized by that at 2 :igm, 1 = 0, is shown for three values of Aobs 
as indicated. Here we set 2 igm,u = 6.295. 
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Fig. 8. The same as Fig. 2, but showing some different models. The three solid curves are the model absorption by the DLA at 
^DLA = 6.295, but with different values of log Ahi = 20.0,20.7, and 21.62, as indicated. The short-dashed curve shows the model 
absorption by the IGM with 2 :igm,u = 6.295 and xhi = 1-0, and the long-dashed curve is the same but the DLA absorption with 
log Ahi = 20 and ^dla = 6.295 is added. 


























